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REFLECTIVE  BUTLER  MATRICES 


1.  INTRODUCTION 

A conventional  Butler  matrix  can  be  formed  by  use  of  symmetric  3-dB  hybrid  couplers 
[1],  Each  matrix  can  be  used  as  a building  block  for  the  next  larger  one.  The  configu- 
rations are  symmetric  about  a longitudinal  axis  which  divides  the  input  and  output  ports 
into  two  halves,  but  except  for  the  two  smallest  ones,  they  are  not  symmetric  about  an 
axis  midway  between  input  and  output  ports.  Therefore,  these  matrices  act  as  lenses  in 
one  direction.  That  is,  the  feed  points  correspond  to  a focal  plane  and  the  outputs  corre- 
spond to  the  region  of  collimation.  However,  because  of  the  lack  of  symmetry,  the  lens 
characteristic  does  not  hold  in  the  opposite  direction;  the  output  ports  do  not  correspond 
to  a focal  plane  with  attendant  collimation  properties  on  the  side  of  the  input  ports.  It 
is  evident  that  a matrix  symmetric  about  an  axis  midway  between  the  input  and  output 
ports  is  desirable.  Such  a matrix  might  be  cut  in  that  plane  of  symmetry  so  that  the 
input  and  output  ports  would  be  identical.  In  this  manner,  the  half-matrix  would  corre- 
spond to  a reflection-type  system  in  which  the  feed  positions  are  in  the  plane  of  the 
aperture.  Thus,  the  disadvantage  incurred  in  achieving  the  saving  of  one-half  of  the  com- 
ponents would  be  the  necessity  for  switching  the  output  between  the  feeds  and  the  radi- 
ators or  possibly  avoiding  feed  blockage  by  the  use  of  circulators. 

In  this  report,  the  procedure  for  modifying  a conventional  Butler  matrix  to  achieve 
a reflective  network  is  presented. 

I 

2.  NETWORK  AND  RADIATION  PATTERN  CHARACTERISTICS  ! 

OF  A BUTLER-MATRIX-FED  ARRAY 

For  the  reader  who  is  unfamiliar  with  the  operation  of  a Butler  matrix,  a brief 
review  is  given  here.  A Butler  matrix  is  a 2N-port  network,  where  N = 2^  (p  is  an  integer). 

All  ports  are  matched,  and  the  N ports  on  one  side  are  mutually  isolated,  as  are  the  N 
ports  on  the  other  side.  The  power  transfer  coefficient  between  any  port  on  one  side  and 
any  port  on  the  other  side  is  1/N.  One  set  of  N ports  is  termed  the  inputs,  and  the  other 
set  is  termed  the  outputs.  If  power  is  fed  into  any  of  the  input  ports,  it  is  split  uniformly 
among  the  N output  ports,  without  loss.  For  each  input  port  used,  there  will  be  a partic- 
ular phase  distribution  among  the  output  ports.  For  appropriate  ordering  of  the  output 
ports,  all  of  the  phase  distributions  are  linear;  that  is,  if  the  output  ports  are  numbered 
n = 1,  2,...,  N,  the  phase  difference  between  ports  n and  n - 1 is  constant  for  all  n.  This 
constant  is  different  for  each  input  port.  If  the  input  ports  are  numbered  m = 1,  2,...,  N,  the 
phase  difference  can  be  expressed 

Mm  = <^o  + 
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where  0o  is  a selectable  constant,  which  is  fixed  for  all  m. 

The  transfer  phase  from  port  m to  port  n can  be  expressed  as 

* 2vm/N), 

v'here  is  a selectable  constant  for  each  value  of  m. 

The  value  of  0o  is  generally  determined  by  the  desired  application  for  the  network. 
For  example,  matrix-fed  circular  arrays  require  cyclic  output  phase  distributions  for  which 
<I>Q  = 0.  The  networks  presented  in  this  report  have  0q  = rr/Af  as  a result  of  the  symmetry 
imposed  upon  them.  It  is  noted  that  can  be  altered  to  any  desired  value  for  any  given 
network  simply  by  adding  an  appropriate  set  of  linearly  increasing  phase  shifts  to  the 
output  ports. 

If  the  output  ports  are  connected  to  a linear  antenna  array,  the  uniform  amplitude 
and  phase  distribution  generated  by  any  input  will  produce  a directive  radiation  pattern 
of  the  form 

8inAr(|i-Po) 

N sm(n-Ho)  ’ 

where 

ju  = 2ird  sin  0/X, 
with 

d = interelement  spacing, 

X = wavelength, 

6 = angle  relative  to  normal  to  the  array, 

and  where 

Po  = 2itd  sin  00  A = 

Oq  = beam  direction. 

The  resulting  multiport  antenna  system,  consisting  of  an  AT-element  array  connected 
to  a 2N-port  Butler  matrix,  has  N ports,  each  of  which  produces  a beam.  It  can  be  diown 
that  the  multiple  antenna  patterns  form  an  orthogonal  set,  as  do  the  output  distributions 
to  the  array.  A beam  corresponds  to  an  aperture  distribution  and  also  to  a particular 
input  port. 

The  network  with  the  above  characteristics  which  uses  a minimum  number  of 
components  is  a Butler  matrix.  This  network  is  closely  analogous  to  the  fast  Fourier 
transform  [2].  In  the  analysis  presented  here,  the  networks  will  be  generated  using 
3-dB  hybrid  directional  couplers,  although  other  types  of  four-port  hybrid  networks  can 
also  be  used. 
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3.  CONFIGURATION  OF  A BUTLER  MATRIX 


Figure  1(a)  shows  the  configuration  of  a 3-dB  hybrid  directional  coupler,  which  has 
four  ports  and  N = 2.  Inputs  1 and  2 form  two  beams  in  the  directions  (Iq  = ±90°.  A 
network  with  N = 4 can  be  formed  from  the  fundamental  hybrid  coupler  as  shown  in 
Fig.  1(b).  Four  beams  can  be  formed,  and  they  are  pointed  at  Hq  = ±45°,  ±135°.  In 
general,  networks  with  larger  numbers  of  input  ports  can  be  formed  from  networks  with 
fewer  ports.  As  an  example,  we  shall  show  how  a 16-port  Butler  matrix  can  be  formed 
from  4-port  matrices.  First,  we  recognize  that  a four-port  matrix  (N  = 4)  transforms  four 
inputs  into  four  distinct  beams  evenly  distributed  in  p space.  Therefore,  we  initially 
require  four  blocks  (each  block  is  a four-port  matrix)  connected  as  shown  in  Fig.  2(a). 
Outputs  of  block  1 become  outputs  1,  5,  9,  and  13,  and  those  of  block  2 become  outputs 
2,  6,  10,  and  14,  etc.  An  additional  row  of  four  4-port  matrices  is  required  to  combine 
the  inputs  of  these  four  blocks.  The  connection  of  these  matrices  to  the  first-row  blocks 
is  shown  in  Fig.  2(b).  The  first  inputs  of  the  four  blocks  in  the  first  row  are  connected 
to  block  1 of  the  second  row,  and  four  inputs  are  formed:  1,  5,  9,  and  13.  Since  the 
beam-pointing  directions  (or  phase  gradients)  of  the  16-port  matrix  are  different  from 
those  of  the  4-port  matrix,  additional  phase  shift  must  be  provided  to  make  up  this  differ- 
ence. For  example,  input  1 of  the  block  has  a phase  gradient  of -45°,  but  in  a 16-port 
matrix,  the  phase  gradient  of  the  first  input  is -11. 25°.  To  account  for  this  difference, 
phase  shifts  of  33.75°,  67.5°,  and  101.25° are  required  for  the  second,  third,  and  fourth 
output  lines  of  the  first  block  in  the  second  row,  respectively.  The  inputs  of  this  block 
now  form  beams  pointing  in  the  directions  Hq  =-11.25°,  -101.25°,  168.75°,  and  78.75°. 
Similar  additional  phases  are  required  for  blocks  2,  3,  and  4 of  the  second  row.  The 
input  lines  of  the  overall  network  are  then  interwoven  to  be  symmetric  with  the  output 
lines. 
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(b)  Four-port  network 


Fig.  1 —Butler  matrices  for  N - 2 and  Af  = 4 
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In  general,  the  matrix  network  layout  and  connections  can  be  described  by  the 
following  procedures  [3] : 

• The  product  of  the  number  of  matrices  in  each  row  and  the  number  of  ports  in 
each  matrix  block  is  equal  to  the  total  number  of  ports. 


• Let  the  rows  be  numbered  1,  2,.,.,  P.  Let  the  blocks  in  the  ith  row  have  L,  ports. 

The  total  number  of  ports  in  the  matrix  is  AT  = IT  L-  and  the  number  of  blocks  in  the 

P,  1=1  ' 

ith  row  is  = .[1 1,,-,  where  the  prime  indicates  that  the  factor  for  ; = / is  omitted. 


• A variety  of  procedures  can  be  established  for  interconnecting  the  blocks  in 
successive  rows  of  the  overall  matrix,  exactly  analogous  with  the  variety  of  flow  diagrams 
available  for  the  fast  Fourier  transform  [4] . A simple  procedure  is  described  here:  The 
rows  to  be  connected  are  i and  / + 1;  the  block  sizes  are  L,-  and  and  the  numbers 
of  blocks  in  each  row  are  M,-  and  where  L,M,-  = = N.  The  outputs  of  the 

blocks  in  row  i + 1 are  connected  with  the  inputs  of  the  blocks  in  row  i,  which  are  num- 
bered sequentially  1 through  N.  The  outputs  of  the  first  block  in  row  j + 1 are  connected 
to  inputs  1,  Mj+i  + 1,  + 1,  etc.  The  outputs  of  the  second  block  are  connected 

to  inputs  2,  i + 2,  + 2,  etc.  In  general,  the  outputs  of  the  ;th  block  in  row 

I + 1 are  connected  to  inputs;  + where  k = 0 through  i - 1. 


• It  is  now  necessary  to  determine  the  transmission-line  phase  shifts  which  must  be 
added  to  ensure  beam  collimation  for  all  input  ports.  For  the  blocks  in  the  first  row, 
which  have  ports,  the  beam-pointing  directions  when  one  of  these  blocks  is  connected 
to  an  Lj  -element  array  (if  l-iis  even)  aie  Hq  = ± jr/Lj,  ±3n/Lj , etc.,  with  a spacing  in  n of 
2ir/Li.  If  Li  is  odd,  (jlq  = 0,  ±2ir/Li,  tAir/Ly,  etc.*  Since  the  output  ports  of  the  first-row 
blocks  are  fed  to  elements  with  spacings  increased  by  a factor  of  A/j , the  beam  maxima 
are  at  Hq  = +3ir/MjLj,  etc.,  or  Hq  = ±ir/N,  ±3ii/N,  etc.;  it  is  also  seen  that 

the  beams  from  a single  block  exhibit  grating  lobes.  Appropriate  phase  shifts  must  be 
inserted  in  the  interconnecting  lines  so  that  a selected  grating  lobe  is  reinforced  by  con- 
tributions from  all  blocks  in  a row  and  all  other  grating  lobes  are  suppressed.  As  in  the 
case  of  the  interconnecting  line  geometry,  there  are  various  procedures  for  achieving  this 
result.  Throughout  this  report  we  have  chosen  to  generate  an  innermost  beam  from  an 
outermost  matrix  port.  The  general  procedure  for  determining  the  phase  shifts  to  be 
inserted  in  the  output  lines  of  any  given  block,  after  one  selects  a beam  position  for  a 
given  overall  matrix  input  port,  is  to  set  down  the  required  output  phase  distribution  for 
the  overall  matrix.  Then  the  block  outputs  in  any  given  row  can  be  independently  ad- 
justed for  the  very  important  reason  that  the  net  transmission  phases  for  all  outputs  of 
any  block  through  the  remaining  rows  of  the  matrix  are  identical. 


^is  procedure  has  been  programmed  in  FORTRAN  language,  and  the  program  is 
described  and  listed  in  the  Appendix.  Sample  computer  plots  are  shown  in  Fig.  3,  with 
N = 32,  Lj  = 4,  L2  = 2,  L3  = 4,  and  in  Fig.  4,  with  N = 64,  Lj  = - 8.  Since  the 


*Although  N is  a power  of  two  for  the  conventional  Butler  matrix,  we  assume  in  this  discussion  that 
the  factors  of  N,  can  be  arbitrarily  chosen. 
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matrix  is  symmetric  about  a longitudinal  tu.is,  only  the  left  half  of  the  matrix  is  shown. 
Tables  1 and  2 list  the  phase-shift  values  that  must  be  inserted  in  the  outputs  of  the 
indicated  rows  of  the  matrices  shown  in  Figs.  3 and  4,  respectively.  The  phase-shift  values 
are  listed  in  the  tables  from  left  to  right,  corresponding  to  the  output  ports  of  the  indi- 
cated rows,  taken  in  order  from  left  to  right  in  the  figures.  Figiires  3 and  4 plot  the  left 
halves  of  right-left  symmetric  matrices,  and  Tables  1 and  2 list  phase-shift  values  for  full 
matrices.  Note  the  symmetry  of  the  phase-shift  listings. 


OUTPUTS 


I — i - I-  1 'i  1 ‘ I I 1 

INPUTS 

Fig.  3 — Left  half  of  a 32-port  Butler  matrix 
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Table  1 — Phase  Shifts  for  32-port  Matrix, 
Inserted  in  Outputs  of  Indicated  Rows. 


ROW  2 

0.0000 

67.5000 

0.0000 

22.5000 

22.5000 

0.0000 

67.5000 

0.0000 

0.0000 

67.5000 

0.0000 

22.5000 

22.5000 

0.0000 

67.5000 

0.0000 

0.0000 

67.5000 

0.0000 

22.5000 

22.5000 

0.0000 

67.5000 

0.0000 

0.0000 

67.5000 

0.0000 

22.5000 

22.5000 

0.0000 

67.5000 

0.0000 

ROW  3 

0.0000 

39.3750 

78.7500 

118.1250 

0.0000 

28.1250 

56.2500 

84.3750 

0.0000 

16.8750 

33.7500 

50.6250 

0.0000 

5.6250 

11.2500 

16.8750 

16.8750 

11.2500 

5.6250 

0.0000 

50.6250 

33.7500 

16.8750 

0.0000 

84.3750 

56.2500 

28.1250 

0.0000 

118.1250 

78.7500 

39.3750 

0.0000 

Table  2 - 

Phase  Shifts  for 

64-Port  Matrix, 

Inserted 

in  Outputs  of  Second  Row. 

0.0000 

19.6875 

39.3750 

59.0625 

78.7500 

98.4375 

118.1250 

137.8125 

0.0000 

14.0625 

28.1250 

42.1875 

56.2500 

70.3125 

84.3750 

98.4375 

0.0000 

8.4375 

16.8750 

25.3125 

33.7500 

42.1875 

50.6250 

59.0625 

0.0000 

2.8125 

5.6250 

8.4375 

11.2500 

14.0625 

16.8750 

19.6875 

19.6875 

16.8750 

14.0625 

11.2500 

8.4375 

5.6250 

2.8125 

0.0000 

59.0625 

50.6250 

42.1875 

33.7500 

25.3125 

16.8750 

8.4375 

0.0000 

98.4375 

84.3750 

70.3125 

56.2500 

42.1875 

28.1250 

14.0625 

0.0000 

137.8125 

118.1250 

98.4375 

78.7500 

59.0625 

39.3750 

19.6875 

0.0000 



4.  SYMMETRIZING  THE  BUTLER  MATRIX 

The  modification  of  the  matrices  to  a symmetric  reflection  form  is  accomplished  in 
three  phases.  First,  the  network  must  be  rearranged  so  that  the  hybrid  couplers  and  in- 
terconnecting lines  are  symmetric  about  the  midline  of  the  matrix,  without  regard  to  any 
of  the  phase  shifters  in  the  interconnecting  lines.  Second,  a scheme  for  modifying  the 
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phase  shifters  to  achieve  symmetry  about  the  midline  is  required.  Finally,  a method 
whereby  the  matrix  can  actually  be  cut  in  two  is  necessary. 

The  method  of  achieving  topological  symmetry  is  based  on  the  requirement  for  a 
fixed  arrangement  of  .he  output  ports,  determined  by  the  output  phase  distributions, 
with  the  accompanying  “natural”  arrangement  of  the  output  couplers.  This  allows  us  to 
arrange  the  inputs  and  input  rows  of  couplers  in  the  same  fashion,  as  shown  in  Figure  5(a) 
for  the  8-port  matrix.  The  next  step  is  to  determine  the  locations  for  the  middle  row  of 
hybrids.  This  is  accomplished  by  averaging  the  positions  of  the  four  hybrids  to  which 
each  middle-row  hybrid  is  connected.  The  result,  shown  in  Figure  5(b),  is  two  centered 
hybrids  and  two  off-center.  The  network  is  now  seen  to  be  topologically  symmetric. 

It  is  necessary  to  modify  the  phase  shifts  so  that  the  network  is  completely 
symmetric  about  the  midline.  It  is  first  required  to  find  an  additional  degree  of  freedom 
beyond  what  was  used  in  the  conventional  Butler  matrix;  that  is,  all  phase  shifts  were 
inserted  in  the  interior  of  the  network— but  there  is  no  reason  why  additional  phase 
shifts  cannot  be  inserted  in  both  outputs  of  any  of  the  output  couplers  and  correspond- 
ingly in  both  of  the  input  transmission  lines  to  any  of  the  input  couplers. 


Fig.  5a  —Symmetric  arrangement  of  input  and  output  couplers  and  ports 
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Fig.  5b  Location  of  middle  couplers  and  symmetric  layout  of  eight-port  matrix 


With  this  added  degree  of  freedom,  it  is  now  possible  simply  to  average  the  symmet- 
rically located  phase  shifts  in  the  interior  of  the  matrix  of  Figure  5(b).  When  a pair  of 
phase  shifts  are  averaged,  it  is  necessary  to  determine  the  value  of  the  compensating  phase 
shift  to  be  added  to  the  outputs  and  inputs  of  the  matrix.  Since  the  relative  output  phase 
distribution  is  fixed,  one  needs  only  to  calculate  the  phase  shift  of  any  one  output  line 
from  a block,  and  the  remaining  output  line(s)  from  that  block  will  be  automatically 
correct  with  the  same  phase  shift.  It  is  also  noted  that  adding  a phase  shift  to  an  input 
line  of  a block  does  not  alter  the  relative  phase  gradient  produced  by  that  block.  In  the 
case  of  the  eight-port  network,  this  compensating  phase  shift  is  +11.25°  inserted  in  the 
outputs  of  the  second  and  third  blocks  of  the  first  row.  The  resulting  completely  symmet- 
ric network  is  shown  in  Fig.  5(c). 

The  objective  with  larger  matrices  is  to  generate  them  by  some  procedure  that  does 
not  increase  in  complexity  as  the  matrix  increases  in  size.  For  example,  a conventional 
Butler  matrix  with  N ports  can  be  formed  from  two  of  the  next  smaller  matrices  with 
N/2  ports  plus  a row  of  hybrid  couplers.  It  is  also  possible  to  devise  building-block  tech- 
niques for  symmetric  matrices.  Three  arrangements  for  a 16-port  matrix  are  shown  in 
Fig.  6.  Figure  6(a),  with  all  Lf  = 2,  is  the  symmetric  layout  for  four  rows  of  hybrid 
couplers;  not  only  was  the  layout  difficult  to  deduce,  but  there  are  two  or  three  degrees 
of  freedom  in  the  phase-shift  adjustment.  In  Fig.  6(b),  two  rows  of  four-port  networks 
are  used.  The  matrix  of  Fig.  6(c)  has  a middle  row  of  four -port  networks  and  top  and 
bottom  rows  of  hybrid  couplers. 
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In  general,  any  Butler  matrix  can  be  realized  by  either  a two-row  configuration 
(Fig.  6(b))  or  a three-row  configuration  (Fig.  6(c)).  The  two-row  method  is  applicable 
to  networks  with  2^"  ports.  For  networks  with  2^"*^  ports,  three  rows  are  required,  and 
the  first  and  third  rows  can  use  couplers:  Li  = = 2.  For  a two-row  configuration,  the 

topology  is  automatically  symmetric.  For  both  the  two-row  and  three-row  cases,  one  set 
of  compensating  phase  shifts  can  be  added  to  establish  symmetry.  In  the  two-row  case,  it 
is  inserted  in  the  lines  connecting  the  rows  of  blocks  as  shown  in  Fig.  6(b).  In  the  three- 
row  case,  the  averaging  procedure  is  used. 

5,  CUTTING  THE  MATRIX  TO  MAKE  IT  REFLECTIVE 

Now  that  a method  for  generating  symmetric  Butler  matrices  has  been  found,  it  is 
desirable  to  devise  a method  for  cutting  them  in  two  along  the  midline  in  order  to  reduce 
the  number  of  components  by  a factor  of  two.  At  first  glance,  one  might  conclude  that 
such  a cut  is  impossible,  because,  in  general,  one  part  of  the  resulting  half -matrix  is  inevi- 
tably isolated  from  another  part  if  the  cut  is  a simple  break  in  transmission  lines  or  com- 
ponents. However,  let  us  first  examine  what  can  be  done  with  the  simplest  matrices  and 
then  perhaps  attempt  to  extend  the  procedure  to  the  more  complex  ones.  The  simplest 
matrix  is,  of  course,  the  3-dB  hybrid  coupler  itself.  The  midline  of  symmetry  passes 
through  the  center  of  the  coupler,  and  it  is  here  that  the  cut  would  be  made.  The  question 
is  what  constitutes  half  of  a 3-dB  coupler?  It  is  seen  in  Fig.  7 that  a 3-dB  directional  cou- 
pler can  be  replaced  by  two  cascaded  8.3-dB  directional  couplers,  with  the  result  that  the 
midline  of  symmetry  now  crosses  the  two  transmission  lines  joining  the  8.3-dB  couplers. 
Therefore,  the  half-matrix  for  the  simplest  case  is  simply  an  8.3-dB  coupler  with  outputs 
open-circuited. 


924  .1.383 

1 I 

i j S-S-dB 
\y  DIRECTIONAL 
/\  COUPLER 

It 


(a)  Operation  of  3 dB 
directional  coupler 


/r 


E 


^DUNE 

OFSYMMETRY 


- OPEN  CIRCUIT 


2 1 


,E 


8.3^B  COUPLER 


(b)  3-dB  coupler  made 
from  two  8.3  dB 
couplers 


(c)  Network  equivalent  to 
half  of  3-dB  coupler 


Fig.  7 —3-dB  coupler  realized  as  symmetric  cascade 
of  two  8.3-dB  couplers 
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If  we  consider  next  the  four-port  Butler  matrix,  shown  in  Fig.  8(a),  the  midline  of 
symmetry  for  this  case  passes  through  the  interconnecting  transmission  lines  between 
directional  couplers.  It  is  seen  that  if  these  lines  are  simply  cut,  the  result  is  two  isolated 
3-dB  couplers  with  open-circuited  outputs.  It  is  found  that  the  only  alternative  is  to  join 
the  two  transmission  lines  that  cross  at  the  midline,  as  shown  in  Fig.  8(b). 


(a)  Symmetric  matrix  with  (b)  Haif-matrix  resuiting  from 

desired  cut  indicated  joining  iines  at  cut 


Fig.  8 — Half  four-port  matrix  illustrating  method  of  making  cut 


The  general  rule  that  emerges  from  these  simple  examples  is  that  transmission  lines 
and  couplers  are  cut  on  the  midline  of  symmetry  of  the  matrix  but  any  lines  or  com- 
ponents which  are  crossed  or  superimposed  on  the  midline  are  not  open-circuited  but  are 
joined  together.  Further  results  of  the  application  of  this  rule  are  shown  in  Fig.  9 for 
N=  8,  in  Fig.  10  for  N = 32  with  Lj  = Lg  = 2,  Lg  = 8,  and  in  Fig.  11  for  AT  = 64  with 
I<1  = Lg  = 8.  Tables  3 and  4 list  the  phase  shifts  for  the  32-port  and  64-port  reflective 
matrices,  respectively.  The  phase-shift  values  are  listed  in  the  tables  from  left  to  right, 
corresponding  to  the  output  ports  of  the  indicated  rows,  taken  in  order  from  left  to  right 
in  the  figures.  Figures  10  and  11  plot  the  left  halves  of  right-left  symmetric  configurations 
and  Tables  3 and  4 list  phase-shift  values  for  both  halves,  with  resultant  symmetry  in  the 
listings.  Since  7/8  of  the  inputs  to  row  1 of  Fig.  11  are  interconnected,  the  values  listed 
in  Table  4 reflect  those  interconnections.  Thus,  the  phase  shifts  for  lines  2 and  9 are  the 
same,  as  are  3 and  17,  4 and  25,  etc. 

It  is  possible  to  verify  the  operation  of  these  networks  by  tracking  signal  paths  from 
any  input  port  through  the  network  to  all  input  ports.  A result  of  this  reflective  opera- 
tion is  that  1/N  ot  the  input  power  returns  to  the  input  port,  thereby  mismatching  it. 

Thus,  the  input  voltage  standing-wave  ratio  (VSWR)  of  a reflective  matrix  is  given  by 

0-V3V+  1 

■ 
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Table  3 — Phase  Shifts  for  Reflective  32-Port  Matrix. 


ROW  1 

0.0000 

0.0000 

0.0000 

0.0000 

19.6875 

19.6875 

19.6875 

19.6875 

28.1250 

28.1250 

28.1250 

28.1250 

25.3125 

25.3125 

25.3125 

25.3125 

25.3125 

25.3125 

25.3125 

25.3125 

28.1250 

28.1250 

28.1250 

28.1250 

19.6875 

19.6875 

19.6875 

19.6875 

0.0000 

0.0000 

0.0000 

0.0000 

ROW  2 

0.0000 

42.1875 

19.6875 

16.8750 

50.6250 

2.8125 

92.8125 

0.0000 

0.0000 

59.0625 

14.0625 

28.1250 

39.3750 

8.4375 

75.9375 

0.0000 

0.0000 

75.9375 

8.4375 

39.3750 

28.1250 

14.0625 

59.0625 

0.0000 

0.0000 

92.8125 

2.8125 

50.6250 

16.8750 

19.6875 

42.1875 

0.0000 

Table  4 — Phase  Shifts  for  Reflective  64-Port  Matrix 
Inserted  in  Inputs  of  Row  1. 


0.0000 

19.6875 

39.3750 

59.0625 

78.7500 

98.4375 

118.1250 

137.8125 

19.6875 

33.7500 

47.8125 

61.8750 

75.9375 

90.0000 

104.0625 

118.1250 

39.3750 

47.8125 

56.2500 

64.6875 

73.1250 

81.5625 

90.0000 

98.4375 

59.0625 

61.8750 

64.6875 

67.5000 

70.3125 

73.1250 

75.9375 

78.7500 

78.7500 

75.9375 

73.1260 

70.3125 

67.5000 

64.6875 

61.8750 

59.0625 

98.4375 

90.0000 

81.5625 

73.1250 

64.6875 

56.2500 

47.8125 

39.3750 

118.1250 

104.0625 

90.0000 

75.9375 

61.8750 

47.8125 

33.7500 

19.6875 

137.8125 

118.1250 

98.4375 

78.7500 

59.0625 

39.3750 

19.6875 

0.0000 

6.  CONCLUSIONS 

In  this  report,  the  generation  of  a Butler  matrix  from  3-dB  couplers  or  larger  building 
blocks  has  been  discussed.  Such  a matrix  has  been  made  symmetric  by  first  rearranging  the 
blocks  and  then  inserting  appropriate  compensating  phase  shifts.  Finally,  a reflective  con- 
figuration has  been  realized  by  applying  rules  for  cutting  along  the  midline  of  symmetry. 

A program  for  generating  and  plotting  these  matrices  has  been  written. 
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Appendix 


A listing  of  the  computer  program  used  in  this  report  to  plot  both  a conventional 
Butler  matrix  and  a reflective  matrix  is  included  in  this  appendix.  If  the  number  of  ports, 
the  size  of  the  basic  building  blocks  and  the  number  of  rows  are  specified,  this  program  plots 
the  interconnections  and  the  relative  phase  shifts  between  building  blocks.  As  an  example, 
a 16-port  Butler  matrix  network  plotted  by  this  program  is  shown  in  Fig.  6.  To  use  this 
program,  two  data  cards  are  required  to  enter  the  necessary  information.  The  first  data  card 
enters  the  following  fixed-point  (15  format)  variables: 

• NPT  — Number  of  ports  of  the  Butler  matrix  network  to  be  plotted. 

• NROW  — Number  of  rows  of  this  network. 

• LL  — A two-digit  number  to  control  what  type  of  Butler  matrix  to  be  plotted. 

The  first  digit  conveys  the  following  instruction: 

1 — Plot  a conventional  Butler  matrix  network; 

2 — Plot  a symmetric  Butler  matrix  network  (symmetric  about  an  axis 

midway  between  input  and  output  ports); 

3 — Plot  a reflective  Butler  matrix  network; 

0 — Plot  both  a conventional  and  a reflective  Butler  matrix. 

The  second  digit  specifies  how  the  plot  is  presented: 

1 — Plot  the  matrix  network  without  phase  shifts  (phase  shifts  will  be 

printed); 

2 — Plot  only  half  of  the  array  (left  half); 

0 — Plot  the  full  network  with  phase  shifts. 

For  example,  if  one  sets  LL  = 23,  since  the  “units”  digit  is  3,  a reflective 
Butler  matrix  will  be  plotted.  The  “tens”  digit,  2,  specifies  that  only  the  left  half  of  this 
matrix  network  will  be  plotted  (without  phase  shifts). 

• XMS  — A floating-point  number  (use  F10.6  format)  to  specify  the  size  of  the  plot 

(in  the  x direction).  The  y direction  is  limited  by  the  computer  plotter  to  be 
19  cm  (7.5  in.).  If  XMS  is  set  to  be  zero,  the  program  automatically  sets  1.9  cm 
(0.75  in.)  between  each  two  input  ports  and  plots  phase  shifts  between  rows 
(ignoring  the  conditions  set  by  LL).  Thus,  a 16-port  network  will  be  plotted 
at  a size  of  19  by  30.5  cm  (7.5  in.  by  12  in.). 

The  second  data  card  specifies  the  number  of  ports  in  each  basic  building  block  in  each 
row.  For  example,  to  plot  the  network  shown  in  Figure  6(c),  three  numbers  are  required  to 
be  entered  on  the  second  data  card:  2,  4,  2.  Thus,  eight  2-port  blocks  are  used  in  the  first 
and  third  rows  of  this  network,  and  four  4-port  blocks  are  used  in  the  second  row.  This  card 
also  uses  15  fixed-point  format. 
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PROARAm  atmtfft 

COMMON/2/PLTAftRAV (2S*) 

OTMITmSTON  MC(9«5121  *PHA(5«SI2) 
dimension  NRP(t>«  >«NbK(6«  ) 

CALL  PLDTS(PLTARRAr#?5A.l8) 

XMaO. 

LLLMlt  PLOT  conventional  FJLL  MATRIX 
lllb?«  Plot  stmmtricaL  toll  matrix 
lllb3,  mlot  half  Matrix 
LLL»®t  Plot  all 

IF  ll»rt,io  plot  Matrix  without  angle 
IF  ll.gt,  gn  plot  half 

IF  LL.RT,  ^0  plot  halfmatrix  (IN  V OIRFCTION)  WITHOUT  ANGLE 
IF  ANV»LtS  IS  TO  9E  PLOTTED  SET  XM  FOU  L TO  ZERO 
read  lOU»NTP.NROW»LL«XMS 
00  format  (-^ISiFIO. 6) 
lF(F0FiR0>«.3 
3 RFaD  1 'If (NRP(I) tI»I»NHOW) 

101  F()RMaT{16IS) 

LLL»MOOiLl fl") 

Lhc«LL/T0 

IF(XMS.LF,0,)XMS«XM 

X«»XMS 

PRINT  110, NTP 

Ho  FORMAT!//,  1ox,»TOTaI  NUMBER  OF  PORTS  IS*, IS,//) 

Print  hi.nrow 

ni  FORMAT'//, TOX, •NUMBER  OF  RWOS  IS*,IS»//» 

PRINT  1 12, (NBRCI) ,I«l ,NR0W) 

112  format (^Ox, •number  of  PORTS  IN  EACH  dLOCK  OF  FaCH  ROW*,//,10x,5l5) 

IF(KC.«T.'')CALL  plot  (XM,<:,,o,,.B) 

*<C«I<C*1 

NRlaNROW,) 

CALL  NT<»K(NTP,NHI,  NBP,NHK,MC,»HA> 

LPsO 

Call  prtoi)T(ntp,nri ,mc,pha,lp) 

IF(LLL  ,r,T.I)r,0  TO  * 

NHaLHC*10 

Call  ntwkri  t<ntp»nrt , nbp,nrk,nh  ,mc,pha,xm) 
ifilll.'je.dgo  To  1 

Call  PlOT(XM*2.,0,,-3) 

A Call  hlFMTX (NTPfNRl ,NHP,NbX,MC,PHa) 

LPal 

Call  PRTOilT(NTP,NRl,MC«PHA,LP) 
lF(LLL.eT,?)RO  TO  S 
nh8LHC*io 

Call  ntwkplt(nip,nri , nbp,nbk,nh  ,mc,pma,km) 

IF(LLL.'»F.2)fiO  TO  1 
Call  PLOT(XM*?.,0,,-3) 

S NHaLHC*10*) 

Call  ntwkplT(ntp,nri , nrp,nrk,nh  ,mc,pra,xm) 

GO  TO  1 

? Call  stop  p(  ot 

E-'O 
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SURQOUTI»JF  NT«(K(NTP«NHl«  Nap«NRK«MC«PHft) 

C«»#»«ThIS  SllSROUTINt  FINOS  THE  CONNECTION  OF  « Bl)TLE«  MATRIX  OP  FFT 
C GIVEN  THE  NUMBER  OF  RONS  ANP  THF  NUMBER  OF  PORTS  IN  EaCH  BLOCK  IN 
C Each  RON 

C«»«*«COMPILEU  RY  J.  K,  HSIAO 

C*»«»»FIRST  version  is  COMPILED  ON  MAY  1»l976 

c*»«*«NTP,  numrfr  Of  total  input  ports  or  samples 

C*****NROW«  number  of  robs  REOUiRFO  TO  PERFORM  THE  TRANSFORMATION 

C*»***NHP«  an  array  stores  the  number  of  ports  In  Each  block  at  EaCM 
c row.  Each  rloc'  in  a HON  has  the  same  number  oe  ports 
C***»«NHKt  An  array  stores  THt  NUMBER  OF  BLOCKS  IN  E«CH  RON, 

C*»***MCi  A TNO  OImENSIONaL  array  STORES  THE  CONNCTIONS  OF  THF  NFTNORK, 

C first  INOFX  of  the  array  represents  the  number  of  current  RON.  IHF 

C location  of  the  second  index  represents  the  physical  location  of 

C the  previous  RON  NHTLE  Trtt  contents  of  it  is  the  connection  to  Hf 

C CURRENT  RON 

dimension  MC(NHIiNTP) tPHA<NRl*NTR) 
dimension  NFTS(5  2) 

DIMENSION  NBP(6^  )tNBK<64  ) 

c computes  The  numheh  of  POH15  IN  Each  block 

NROWaNPl.) 

PIatI,UISR?AS3«> 

PI?»PI*^. 

NrP(NRI)«I 
NTP2«nTP/P 
DO  lO  TBitNRi 

10  N^K(T)«nTP/NbP(I) 

C***p  NETS  aRWAY  STORES  Tme  LOCATION  OF  IHE  SAMPLES  IN  FaCH  BEaM(OR 

c Frequency  sample),  the  structure  ts  characterized  by  tno  numbers, 

C NTS.NUMBEB  of  time  SAMPLESIOR  input  PORTS)  ANO  NFS,  number  of 

C FREQUENCY  SAMPLESIOR  NUMBER  OF  BEAMS),  FOP  IXAMPle,  nFTS((B-1)* 

C NTS*1)  Is  THF  physical  LOCATION  OF  Thk  FIRST  TIMfc  SaMPLF  TN  THE 

C third  FRFOIIFNCT  SB0UP(  or  of  the  THIRO  BEAM), THIS  IS  OtPRFSFNTtO 

C BY  LMC 

r 

C SET  THF  INITIAL  NFTS  ARBAY 

OJ  11  lal.NTP 

11  NrTS(I)»I 

C*P#*  NTSl  IS  The  previous  values  of  THF  number  of  TIME  SAMPLES (OR  INPUT 

C PORTS) 

C«»**  NTS2  IS  The  CURRENT  VALUE 

C**#*  NFSl  IS  THE  PREVIOUS  VALUE  OF  THF  NUMhER  OF  FREQUENCY  SaMPLFSOR 
C B£  ) 

C»*«#*NFSp' IS  The  current  VALUE 

c 

c 

C SFT  THF  INIAI-  values  OF  NTS  ANO  NFS 

NISlaNTP 
NFS) ■! 

DO  20  la) ,nRI 

C MM  THE  •'•umber  of  blocks  OF  THE  CURRENT  ROW 

C NN,  The  number  of  ports  in  Each  block  of  the  current  row 

MMaNBK (I) 

NNaNBP(l) 

C SET  NTs*?  aNO  NFS2 

NTSPaNTSl /NN 


J 
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NrS2»NTP/NTS2 

c****  THE  ACTUAL  REQUlRtT  PHASE  OPADIE^^T  BETWEEN  SUCCESSIVE  ELEMENT  fOR 
C The  first  ream  IS 
PABaPI/NPS? 

€*••  available  phase  sraOIEnt  for  the  first  beam  in  Each  block  is 

PS6-PI/NN 

KKaO 

C«««  THE  PH*SE  ANOLE  is  STORED  ACCORDING  TO  THE  INDEX  OF  THE  LOWER  ROW 
DO  30  jsI.mM 
MODJbMOOI JtNFSl) 

IF(MOOJ,EO,0)K(OoJ"NFS1 
JJ»< J«1 ) /NES’ 

PA66«PA6*(M00J*2«n 
PASGDsPSG-PAGG 
DO  30  Ksl.NN 
KlaK.l 

kk>kk«i 

LMC»(M00J-1 »*NTSl* (K-11*NTS24JJ 
MCLOC»nFTS(LMC) 

Mr  n ♦HfLoc) ■KK 
IF<KK,LF.NTP?)60  TO  31 
KKTaNlp-KK*l 
PMA(I«KK)aPHA(I«KKII 

60  TO  30 

31  lElPASRD.GT.O.ieO  TO  32 

PHA(I,KK)  aARSlPASliOlPlNN-K) 

60  TO  30 

3?  PHA(I«KK)aPASGO*Kl 
30  continue 

C recording  the  RREUUENCT  sample  Or  beam  position  into  NftS  array 
NTSiaNTS? 

NESlaNES? 

KKaO 

C MNS  TS  THF  number  of  blocks  within  each  group  of  EREOUENCV  samples 
MNSaMH/NTSi 
00  40  jai ,NFSl 
JmOD>MOO( JtMNS) 

IF (jMOn*EO.O) JMOO"MNS 
JJa(J-n/MNS*l 
00  AW  K»i»NTSl 
KKaKK*l 

Aft  nets  «AK) ■ IK-T ) PNESl ♦ ( JMOO-1 ) •NN* JJ 
20  CONTiNUt 
RETURN 
ENO 
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SUBROUTINE  NTMKPLT(NTPtNRl*  NBP«NRKtLH,HCtPH«,XM) 
C NHal,  hale  matrix (IN  X) 

C NHaO,  full  matrix 

C LPTaO  MITM  ANGLES 

C *-PT"1»  mithOut  angCfs 

c lpt"2#  without  angles  ano  half  (IN  Y> 

dimension  MC(NRi*NTP) tPHA(NRI«NTP) 
COMM0N/2/PlTaRRaV(2SA» 

C0MM0N/3/MCT(5l2) 
dimension  NRP(6A  )*NBK(6A) 

DIMENSION  XCR> •XCH2(S1?) 
dimension  LAP(512) 

NHaMOO(LHtlO) 

LPTaLH/lO 

PI»3,Ul5«»?6516 

RACaien./Pi 

NHOMaNRl'l 

Ml-ET».03Sa?, 

Ml ET«A.»HLET/7. 

YMaT.S 

IE (XM.GT.n.)GO  TO  7 
XMa,7SANTP 

7 lE(XM,LT.10,)XMalO, 

XH2aXM/^ 

YHHa.n 

YSPACEaYM/<NROWAl ) 

XSlaXM/(NTP«1 ) 

NPPaNRnW*? 

NNaNRp* 

YHaYM*1 • 

C IE  NH^It  plot  half  THE  MaTRIX(X-SYMMETBY) 

LLL«P 

IF(NH,NE,i)60  TO  5 

c determine  if  the  nRow  is  eyfn  or  odd 

LNHa ( NROW* ) ) /2-NROE/? 

C LNHal.oDOi  LNHaO,  gvEN 

NNa(NH0W*1)/2*l 
LLLaNH*UNH 
5 00  10  T”i«NR2 

lE(I.6T*NNtGD  TO  6 

C FIX  the  PORT  LUCaTION  IN  EaCH  ROW 
lEd.oT.DOO  TO  1 
DO  1 1 J«1,nTP 
n XrMl (j)«(j.l)*XSl 

YllaYR 
60  TO  10 

1 IF(I.LT.NRr)60  to  ? 

DO  12  jal.wTP 

12  XCM2  (j)a(j.i)*XSl 

Y22aYR 
GO  TO  Y 

2 NNBPaNBP(I-l) 

NNgKaNRKd.n 

XRKa,e*XM/NNAK 

XHKSa.flaXRK 

XLSaXgKS/MNRP 
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X1b,1*XM«.1«XBKS 
Yl«YH  ♦.SbyBH 
V2»YR  -,5«V8H 
Y22-Y1 
KK«o 

00  13  JaltNNRK 

IF(LBT.6E.2t«N0.xl.BT.XM2)80  TO  9 
CALL  PlOT<X1.Y2iJ) 

X2BXUxdKS 
Call  plot(xi«yi«2) 

Call  PL0T(x2.vlt2> 

CALL  PL0T(X2«Y2t2) 
lE(LLL.LT,l,OR.I,NE,NN)eo  TO  B 
CALL  0aSH(X1 •Y2tX2tY2) 

GO  TO  0 

CALL  PLOT(XltY2i2) 

XbXUxhkS/P. 

YbYR 

XPb,5»xLS*X1 
00  u K»1.NNBP 
KKbKK«1 
XCM?(^K)BXP 
U XPbXP*XLS 
13  XIsXIaxSK 

c draw  network  connections 

T 00  2o  J«1 ,nTP 

XIbXcmI ( j» 

LMCbMC(T-1 ,J) 

X2bXcM?(LMC) 

XIbXI 

YIlBYll 

1E(LPT,LT,2.0R.XI,LT,XM2)GO  TO  18 

IF <X2.RT.XM2)Gy  TO  PO 

XTbXM 

YIlBYEP(X?,Y22tXl»ri 1 tXM2) ♦V22 
Ig  XTbx2 

Y2TBY2? 

IF(LPT.LT,?.OR.X2.LT.XM2)GO  TO  17 
XTbXM? 

Y?TBYli-YrP(Xl,Yll,x2tY22tXM2) 

IT  YTIbO. 

IF(1,E0,2)yT1b,1 

IF (LPT,GE,?,AN0.X1.GT.XM2) YTIbo, 

CALL  PLOTlXItYIUYTl  ,3) 

IF(I,GT,2)G0  TO 

IF  (lPT,6F,?.ANi;.X1.GT.XM2>G0  TO  A 
Call  PL0T(XI*Y11i2) 

A CALL  PL0TlXT»Y?Tf2> 

IF (I.NF.NRPIGO  to  ‘9 

if(x^.re,xm2,ani),lpt,eq.2)so  to  I9 
IF(I,to.NR?,ANU.X«:.LT.XM2)CALL  PLOT  < X2  , Y22-.  1 ,2) 
CALL  PL0T(x2»Y^2.,1,2) 

>9  IF(LPT,GT,0)(50  to  20 

IF(LPT.GT.O)60  to  20 
ababs<p«a<t-i »LMC»  >*RAC 
lF(A.l.E..001)GO  TO  20 
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lF(I,Nr.Nn?)90  TO  IS 
Y??2«Y?2*vSp8CE-.5»Y8H*2.«HLET 
X»Xl.3.»WLET 
GO  TO  U 

15  Y222«Y22-P,*HLET 
X«x2-J,#WlfT 

16  CALL  NuMrPRCX  tY222,HLET*A*>*.t«HF8.2) 

20  CONTINUE 

C RESET  XCMi  array 

00  21  J»l#NTP 

21  XCMl (J)»Xcm2( J) 

YIUy^ 

10  YR.YR-YSPaCE 

IFtLRl .GE.plcALL  WASH ( XMp, VM*1 . 1 t XMg, vpp-, i ) 

RETURN 

6 Y»YR*YSPACF/7. 

XHMaXM 

IF (LPT,6F,?»XMM«XM? 

CALL  0ASH(0.,YtXNH«V) 

00  31  .I■1,^|TP 

31  LaP(J)>0 
ZbY-^.aHLFT 
NTPK«nTP 

IF(LPT.6F,?)NTPK«N1p/2 
00  30  J«1.NTPK 
1F(LAP(J) .OT.OIGO  to  30 
XUXcMi  (J) 

JJ»MCT(J> 

LAP<JJ)»i 
X?»XCM1 (JJ) 

XbXI 

XT»X2 

YT«Y2 

IF(LPT.LT.?.0R.X2.LF,XM2)G0  TO  32 
XTaXMi; 

YT»YEP(X.Y»XP»Y2»XM2>»Y 

32  Call  PLOT(xltYP»J) 

CALL  PL0T(X  .Y»2) 

CALL  PLOTtXTfYTt2) 

IF(LPT,6T,0)00  10  30 
IF<LNh.6T.O)60  to  30 
LMC>HC(NN,j) 

AbABS (PHA (NNtLMC) ) •RAC 
XbX-3.bWLFT 

CALL  NUNRCR(X»ZtHLET*Ai0.f AHF6t2) 

30  CONTINUE 

IF(LPT.6E,2)CALL  OASH  (X^Zi  yM^l  , 1 .x*<2f  Y) 

RETURN 

CNO 
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SUSnOUTINf  D«SH(Xl»YltX2tV2) 

C this  subroutine  ORXhS  X DASHED  LINE  FROM  A POINT  XI*  Vl  TO  A POINT 

C X?,  V2 

X12«X2-XI 

Y12«Y?-T1 

0«S0RT(A1P**?*Y12**?» 

SO«*T 

2 NN»0/S0 

OObO>NN«SO 

IFlABSrOO) .LE.. 001)00  TO  I 

SOaSO*nU/MM 

60  TO  ? 

1 XINCbS0*X12/0 

YlNC«So*YlP/r» 

X,Xl 

Y«Y1 

DO  10  Tb1*NN»2 
CALL  PLUTrx,Y,3) 

XBX»XiNC 

y:y*yinc 
call  PL0Trx.Y,2) 

XaX*XlNC 

YaY*YlHC 

10  continue 

RETURN 

end 
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SUBROUT I NE  HLFMTX ( NTP t NR 1 1 NRP « MRK t MC « PH A ) 
dimension  NRP(64)«NRKI64) 
dimension  MC (NRl tNTP) tPHA ( NRl t NTP) 

COMM0N/3/Mcr(5*2) 
dimension  aNOISIZ) tATEMPtolE) 

NNbNR1/2 

tL»(NR)«l)/2-NN 

Lt«i  number  of  rows  is  even 

LL>0  NUNRER  of  ROMS  IS  ODD 

CALL  PHASl)M(NRl*N>P*NBP*MC»PHAf AN6) 

DO  lo  TalfNTP 
Ilal 

DO  11  JaltNRl 
KKaMc(J*IT) 
lF(J,EO.NN)KKPaKK 
1I>Kk 

11  CONTINUE 

fiND  thE  joint  Point  then  store  in  mct  array 

00  12  J»1.NN 

FKS«MC(J#KM» 

12  KKaKKS 
MCT(KNP)at(KS 

averagp  The  phase  angles  for  symmetrical  matrix 

00  13  lal.NN 

JJaNRl.j*l 

lMCaMC(J»t» 

AYOa(PHA(J«IMC)*PHA(JU*II )/?• 

>3  PHA<J»TMc)aPMA( JJtllaAVO 
•O  CONTINUE 

CORRECT  phase  ANbLe  of  THE  HIDOLE  ROa  WHEN  THE  NoMrER  OF  ROWS  JS 
Even 

IFILL.LE.oISO  to  I 
NlaNN*) 

DO  20  Tsi.MTP 
I1wMc(n1#T» 

INaHc(Nl*McT(n) 

ATEMPdliaPHAlN  (lit 

IFIPH  (N1  ,IN)  .l*r.AlFMP(II))ArEMP(in«PHAlNl«IN) 

20  CONTINUE 

CORRECT  THE  PHASE  ANQLE  BY  ADDING  THE  SAME  EXTRA  PHASE  TO  EACH 

PORT  IN  A BLOCS 

NMPaNBPIN]) 

NMbbNRKINI) 

DO  2l  I>I,NMR 
IMga(I.l)aNMP 
AAaO, 

00  2M  jaltNMP 

KKaIHe*J 

A>ATEMP (KK) >PHA (Nl tKK) 

IF(a.GT*AA)AA"A 

22  CONTINUE 
lF(AA*LE.n.)  GO  TO  21 
00  23  jwltNMP 
KKalMB«J 

23  PHA(NltKK)aPHA(Nl*KK)«AA 

21  continue 
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RETURN 

CORRECT  THE  PHASE  ANGLES  FOR  THE  CASE  HHEN  THE  NUMRER  OF  ROWS  IS 
000 

1 CALL  PHASUM(NRltNTP«NBP«MCtPHAtATEMP) 

DO  30  I»IfNTP 
AAaAN6(I)*ATENP(I) 

JJ-MC  (l»T) 

30  PMA(1*JJ)bPHA<NR1*I)sAA 
RETURN 
End 


SUBROUTINE  PHASUM ( NR 1 *NTP « NBP t HC t PHA t AS) 
dimension  NRP(6AI 

dimension  MC(NR1*NTP) »PHA(NRt«NTP) tAS(NTP) 
DIMENSION  LAP(^*&^?)*AIP*^) 

C SET  the  PM*SF  shift  OF  THE  BOTTOM  RO" 

NR0MaNRl»1 
NNaNBP(NROH) 

DO  1 U«I*NN 
LAP<2#J)aj 
I ASlJ)apHA(1*I) 

KKaNN 

00  in  lal.NROM 
IlaNROW.I 
Il>II«t 
NNaNBP(II) 

I^III.LEtO)NN»l 
DO  1?  jal.NTP 
LaP(1  t j)a(.AP(2tJ) 

12  A(J>aAS(J) 

KNa 

00  20  L>lfKK 

LLaLAP(ltL) 

00  21  NaltNTP 
IP(MC(li«N),NE.LL)00  TO  21 
JJaN 

00  TO  ?2 
21  continue 

2?  NMDaMOO(JjtNN) 

IF  <NMO*EO.0)nmO*NN 
00  30  Kal ,NN 
lN0ajJ-NMn*K 
lF(NN.EO«i)INOaJJ 
KNaKNA) 

LAP(2«KN)aTNn 

30  aS(IN  )>A(LLI*PHa(I1 *LL) 

20  CONTINUE 
KK.Kn 

10  continue 

RETURN 

E^n 
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SUdROUT INE  PRT out ( NTR • nr 1 * RC  t PH a * LP ) 

DIMENSION  MC(NRltNrP)«PHAINRl*NTP) 

DIMENSION  PaNOIS  2) 

PIa3,Ui5«P6S36 
RACalRO./pl 
tP(LP.eT.O)  60  TO  I 
PRINT  10* 

102  FORMAT (//tl OX **C0NVENTI0NAU  matrix  CONNECTION  AND  PHASE**//) 

00  TO  F 

1 print  lOT 

lOT  FORMAT(//*10X**STMMETRICAL  MATRIX  PHASE  ANGLES**//) 

00  TO  R 

2 00  10  1*1*NR1 

PRINT  l03,I 

103  FORMAT«//*20X**ROM-**I5*//) 

10  PRINT  lp4*(MCiI*Jl*J*l*NTP) 

F0RMAT(*H  , 10X*1®I5) 

PRINT  ’05 

105  format (//*10X**CONNpcTlON  PHASE  ANeLCS**//) 

3 00  20  I”1*NR1 

DO  21  J«1*NTP 

LMc*J 

IF  < 1 .EO.NRl ) LMC>MC ( 1 * J) 

21  PaN0(J)*PH8(I*LHc)*RaC 
PRINT  l®3,t 

20  PRINT  l06,(P«N0|J)*Jal*NTP) 
lOft  F0RMAT(IH  , 10x*5F10,A) 

Return 

end 


function  VEP(XltYl*X2*Y2*XM) 

XR«(XM-X1)/(X2-XM) 

YEP«XR*aRS(Y1-T2)/(1,*XR) 

return 

end 
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